i?-values in Low Energy e + e Annihilation 



Zhengguo Zhao 

Institute of High Energy Physics of Chinese Acadmy of Sciences 

O ! Beijing 100039,P.R.C. 
O 

: February 7, 2008 

0\ 



Abstract 



' This presentation briefly summarizes the recent measurements of R- values 

in low energy e + e~ annihilation. The new experiments aimed at reducing the 
uncertainties in R- values and performed with the upgraded Beijing Spectrom- 

' eter (BESII) at Beijing Electron Positron Collider (BEPC) in Beijing and with 

CMD-2 and SND at VEEP-2M in Novosibirsk are reviewed and discussed. 

O 
O 
O 

X : 1 Introduction 

p 

According to the quark-parton model, hadrons produced via e + e collisions are 
^ ; characterized by the annihilation of e + e~ pairs into a virtual 7 or Z° boson. In the 

lowest order, one defines the ratio of the rate of hadron production to that for muon 
^ ■ pairs as 

£ : ^^ " ( ff^ ha f om) =3EQ 2 f , (1) 

a(e + e~ — > J j 

where Qf is the fractional charge of the quark, and the factor of 3in front counts 
the three colors for each flavor. <r(e + e~ — > = 47ra 2 /3s is the cross section of 

the pure QED process. The value of R, which counts directly the number of quarks, 
their flavor and colors, is expected to be constant so long as the center-of-mass (cm) 
energy of the annihilated e + e~ does not overlap with resonances or thresholds for 
the production of new quark flavors. One has 

R = 3[(2/3) 2 + (l/3) 2 + (2/3) 2 ] = 2 for u, d, s 
= 2 + 3(2/3) 2 = 10/3 for u,d,s,c, 
= 10/3 + 3(l/3) 2 = 11/3 for u,d,s,c,b. 



These values of R are only based on the leading order process e + e~ — > qq. 
However, one should also include the contributions from diagrams where the quark 
and anti-quark radiate gluons. The higher order QCD corrections to R have been 
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calculated in complete 3rd order perturbation theory ||, and the results can be 
expressed as 

R = 3EQ 2 f [l + (^) + 1.411(^) 2 - 12.8(^) 3 + ...], (2) 

^ 71 7X 7X 

where a s (s) is the strong coupling constant. Precise measurement of R at higher 
energy can be employed to determine ct s (s) according to eq.2, which exhibits a 
QCD correction known to O(af) . In addition, non-perturbative corrections could 
be important at low cm energy, particularly in the resonance region. 

R has been measured by many laboratories in the energy region from hadron 
production threshold to the Z° pole and recently to the energy of W pair pro- 
duction 0. The experimental R- values are in general consistent with theoretical 
predictions, which is an impressive confirmation of the hypothesis of the three color 
degrees of freedom for quarks. The measurements of R in the low energy region were 
performed 15 to 20 years ago in Novosibirsk, Orsay, Frascati, SLAC and Hamburg 
[§. [|, ||, ||, |7], ||, ||. Fig. |l] |jl6| shows the R- values for cm energies up to 10 GeV, in- 
cluding resonances. For cm energies below 5 GeV, the uncertainties in .R-values are 
about 15% on average; and the structure in the charm threshold region is not well 



determined. The DASP group 1X3 1 inferred the existence of narrow resonances at 



4.04 GeV and 4.16 GeV. In addition to the resonance at 3.77 GeV, Mark I data [14 



shows broad enhancements at 4.04, 4.2 and 4.4 GeV. The resonance at 4.4 GeV 



was also observed by PLUTO [15], but the height and width of the resonance were 
reported differently. A new cross section measurement in the charm threshold re- 
gion is needed to clarify the structure, which is important not only for the precision 
determination of a(M§) and the interpretation of the (g — 2)^ measurement of E821 
at BNL but also for the understanding of charmonium itself. 

Between the charm and bottom thresholds, i.e., about 5-10.4 GeV, R was mea- 
sured by Mark I, DASP, PLUTO, Crystal Ball, LENA, CLEO, CUSB, and DESY- 
Heidelberg collaborations. Their systematic normalization uncertainties were about 
5-10%. Above bottom threshold, the measurements were from PEP, PETRA and 
LEP with uncertainties of 2-7%. 

Remarkable progress has been made in the precision test of the Standard Model 
(SM) during the last decade. The electroweak data from LEP is so copious and 
precise that one can make use of the radiative correction effect to test the SM. In 
particular, the indirect determination of m# depends critically on the precision of 
a(M§). Recently, for example, there has been an increasing interest in electroweak 
phenomenology to reduce the uncertainty in a(M§) which seriously limits further 
progress in the determination of the Higgs mass from radiative corrections to the 
SM [TE, T7_, IIS], [BI ^U], ^TJ. The uncertainty in a(M§) arises from the contribution 
of light quarks to the photon vacuum polarization Aa(s) = — n 7 ( s ) & t the Z mass 
scale. They are independent of any particular initial or final states and can be 
absorbed in a(s) = a/[l — Aa(s)], where the fine-structure constant a=l/137.035 
989 5(61). Aa = Aai ep t on + Aa had, °f which the leptonic part is precisely calculated 
analytically according to perturbation theory because free lepton loops are affected 
by small electromagnetic corrections ||12|| . Whereas, the hadronic part Aothad cannot 
be entirely calculated from QCD because of ambiguities in defining the light quark 
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Table 1: R in low energy measured by different laboratories 



Place 


Ring 


Detector 


(GeV) 


Points 


Year 


Beijing 


BEPC 


BESII 


2.0-5.0 


85 


1998-1999 


Novosibirsk 


VEPP-2M 
VEPP-2 


CMD-2 
SND 
01ya,ND CMD 


0.6-1.4 
0.3-1.4 


128 


1997-1999 


SLAC 


Spear 


MarkI 


2.8-7.8 


78 


1982 


Frascati 


Adone 


772, ME A 
Boson,BCF 


1.42-3.09 


31 


1978 


Orsay 


DCI 


M3N,DM1,DM2 


1.35-2.13 


33 


1978 


Hamburg 


Doris 


DASP 
PLUTO 


3.1-5.2 
3.6-4.8 


64 
27 


1979 
1977 
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Figure 1: Experimental R values in the energy region below 10 GeV. The relative 
errors on R in the continuum is given in numbers at the bottom of the figure. 



masses m u and as well as the inherent non-perturbative nature of the problem 
at small energy scale since the free quark loops are strongly modified by strong 
interactions at low energy. An ingenious way to handle this[10] is to relate Acthad 
from the quark loop diagram to R, making use of unitarity and analyticity, 

Aa£ w = _^p, r : ds & + / : d , w) ), (3 ) 

4n z a JAm\ s\s — s) JE^ m s'(s' — s) 

where R(s) = a(e + e~ — > hadrons) / a(e + e~ — > and P is the principal value of 

the integral. 

Much independent work has recently been done to evaluate a(s) at the energy of 
the Z pole. So far, the uncertainty of Aa(s) is dominanted by the i?-values at low 
energy (E cm < 5 GeV) measured with an average uncertainty of ~15%, as indicated 
in Fig. P]. 
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The anomalous magnetic moment of the muon = (g — 2)/2 receives radiative 
contributions that can in principle be sensitive to new degrees of freedom and in- 
teractions. Theoretically, is sensitive to large energy scales and very high order 
radiative corrections. It therefore provides an extremely clean test of electroweak 
theory and may give us hints on possible deviations from the SM. The experimental 
and the theoretical predictions on a M are well reviewed by Dr. Robert Lee in his 
talk given at LP99. 

One can decompose 

a SM = a QED + a ha d + ^ (4) 

where the largest term, the QED contribution a® ED , has been calculated to 0(ct 5 ), 
including the contribution from r vacuum polarization. a™ eak includes the SM ef- 
fects due to virtual W, Z and Higgs particle exchanges. a^ ad denotes the virtual 
hadronic (quark) contribution determined by QCD, part of which corresponds to 
the effects representing the contribution of the running of a(s) from low energy to 
high energy scale. It cannot be calculated from first principles but can be related 
to the experimentally determined R(s) through the expression 



where K(s) is a kernel varying from 0.63 at s — 4m^. to 1.0 at s — oo. 

The hadronic vacuum polarization is the most uncertain of all the SM contri- 
butions to a M , the uncertainty is presently 156 x 10~ n . For several scenarios, it 
has been claimed that "the physics achievement of the effort to measure the cross 
section of e + e~ — ► hadrons that brings down the uncertainty of to 60 x 10 -11 is 
equivalent to that of LEP2 or even the LHC" . 

From equations 3 and 5 one finds that a^ ad is more sensitive to lower energies than 
higher ones. Further measurement in the energy region of 0.5-1.5 GeV from VEPP- 
2M in Novosibirsk and DA$NE in Frascati will contribute to the interpretation of the 



dp measurement at Brookhaven [g(J and the luminosity measurement at CERN |17 
However, their contribution to the precision determination of a(M|) is limited. 
The improved R value from BESII at BEPC in the energy region of 2-5 GeV will 
make the major contribution to evaluate a(M§), and also partly contribute to the 
interpretation of a^. 



2 Recent measurements of R in low energy e + e 

There are two different approaches to the measurement of R. One is to study the 
exclusive hadronic final states, i.e. to measure the production cross section of each 
individual channel a exp (e + e~ — > hadrons)j. The value of R can then be obtained 
by summing over the measured hadron production cross section of all individual 
channels. This method demands that the detector has good particle identification 
and requires the understanding of each channel. It is usually used for cm energies 
below 2 GeV. 
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Another method treats the hadronic final states inclusively. It measures R by 
dealing with all the hadronic events simultaneously and is suitable in an energy 
region where a reliable event generator for hadron production is available. With an 
improved Lund Model ||27||, we may be able to extend this region down to 2 GeV. 



The typical features of hadron production below 5 GeV are: 

• many resonances in this energy region, such as, p, u>, <p, p', u>', </?', cc and charmed 
mesons J/ip, ift(2S), D + D~, D+D~ and r + r~, baryon-antibaryon pair production 

• a small number of final states and low charged multiplicity, usually N c h < 6. 
The experimental challenge here is how to subtract the beam-associated background 
and select N had . 

In the following section, I will first discuss some new measurements done by 
CMD-2 and SND at VEPP-2M in Novosibirsk, which are based on an exclusive 
analysis of the hadron production in the energy region around 0.4-1.4 GeV. Then 
I will concentrate on discussing the R scan done by BESII at BEPC in Beijing in 
the energy region from 2-5 GeV, which measures R values by dealing with hadronic 
final states inclusively. 



2.1 Recent results from VEPP-2M 

VEPP-2M, the e + e~ collider with maximum luminosity of ~ 5 x 10 30 cm _2 s _1 at 
E beam =510 MeV, has been operating since 1974 in the energy region E cm = 0.4 — 1.4 
GeV (p, lu, 0-meson region). SND [^J and CMD-2 j40| are the two detectors 



carrying out experiments at VEPP-2M. Since 1994, VEPP-2M performed a series of 
scans from 0.38 GeV to 1.38 GeV With this data both SND and CMD-2 have 
measured the cross sections for the channels 7r + 7r~, tt + tt~tt°, 7r + 7r~7r + 7r~, 7r + 7r~7r°7r°, 
7r + 7r~7r + 7r~7T° as well as K^Kg and K + K~ . 

• Study of e + e~ — > 7r + 7r~7r + 7r~, 7r + 7r~7r°7r° 

The four pion final states produced via e + e~ annihilation in the E cm = 1 ~ 2 GeV 
energy region dominate and determine the main part of the hadronic contribution 
to a M and the QCD sum rules. Besides, these processes are important sources of 
information for the understanding of hadron spectroscopy, in particular for the study 



of the p-meson radial excitation [31]. These processes were studied at the VEPP- 



2M, DCI and ADONE colliders pi 1331 133, [35], |3§ . The statistical errors of these 



measurements were ~ 5%, and the systematic errors were ~ 15%. There was about 
a 20% discrepancy among the different experiments. 

For the process of e + e~ — > 7r + 7r~7r + 7r~, the systematic error from the recent 
SND results is ~ 7%, mainly coming from the event selection and the luminosity 
determination. The measured total cross sections from CMD-2 for e + e~ — > 27r + 27r~ 
are also illustrated in Fig. 0. Only the statistical errors are shown. The system- 
atic uncertainties are ~ 7%, attributed to the luminosity measurement, the event 
reconstruction and selection and the radiative correction as well. 

For SND, backgrounds to the e + e~ — ► 7r + 7r~7r°7r° channel are mainly from 
e + e~ — > K + K~, QED processes e + e~ — > e + e~e + e~, e + e~77, as well as cosmic 
rays and beam associated background. The systematic error is ~ 7%, of which 
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~ 5% arise out of the variation of the detection efficiency shown from the simula- 
tion of the intermediate states like um°, p°ir°ir and Lorentz-invariant phase space 
simulat ion ( L IP S ) . 

The total cross section for e + e" — ► 7r + 7r~27r° process measured by CMD-2 also is 
plotted in Fig. |||. The error bars indicates only the statistical error. The systematic 
uncertainties mainly come from event reconstruction, radiative corrections and the 
luminosity determination. The overall systematic uncertainty is estimated to be 7%. 
The cross section measured by this experiment is consistent with what was measured 



by OLYA and a recent result from SND ||43fl . However, the cross section from 
all three measurements is apparently lower than that given by ND [53, For 
comparison, the results from Orsay and Frascati above 1.4 GeV are also shown in 
the figure. 

CMD-2 finds that the dominant contribution to the cross section of the process 
e + e" — > 7T + 7r~27r° comes from un° and p^n^n intermediate states, whereas the 
pi^o^o s t a t e i s not observed. The p ± 7i T n° states are saturated completely by the 
ai(1260)7r intermediate state. This is also the dominant contribution to the cross 
section for the process e + e~ — > 2tt + 2tt~ . The theoretical predictions for the differ- 
ential distributions and the total cross sections can be dramatically changed if one 
takes into account the interference of different amplitudes with various intermediate 
states but identical final states. 

The cross section of e + e~ — > Arc can be related to the the four 7r decays of the 
r-lepton through the hypothesis of CVC [|H|]. This has been experimentally tested 



to be valid within an accuracy of 3-5% [41]. The observed ai(1260)7r dominance, if 



it is true, should be taken into account in r decays. 




Figure 2: Energy dependence of the Figure 3: Energy dependence of the 
cross section of e + e~ — > 27r + 27r~ . cross section of e + e~ — > 7i + 7i~2ir°. 



The investigation of e + e — > 7i + n 7r° 
This process was measured by ND at VEPP-2M in the energy region up to 1.1 



GeV |37j • The measured cross section is significantly higher than that predicted by 
the Vector Dominance Model (VDM). However, it is well known that the VDM is 
able to well describe the cross section near the uj and <fi resonances for the processes 
e + e" — > uj, <f) — > pir — > 7r + 7r~7r°. It is therefore necessary to perform new precise 
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measurements in the non-resonance region to investigate the limitation of the VDM 
and determine possible contributions from heavier intermediate states like u;(1120) 
or w(1600). 

SND also measured this channel. The systematic errors from the detection effi- 
ciency, luminosity measurement and the background subtraction are 10%, 5% and 
5% respectively, giving a total systematic error of ~ 12%. The results from the new 
measurement agrees with the old ND data. 

Invariant masses of 7r-meson pairs in the final 3n state were measured to investi- 
gate the intermediate state in the e + e~ — > 7r + 7r~7r° process. The intermediate states 
might be pre and the much less probable ujtt with decay to — > 2n. Comparing the 
mass spectrum of 7r + 7r~ with that of ^n^, one can observe the interference between 
the two intermediate states. The clear peak shown in tt + tt~ mass spectrum of the 
experimental data proves the p — uo interference in 3n final state, and the phase 
measured to be zero agrees with the VDM prediction. Fig. f| plots the cross section 
for the production of 7r + 7r~7r°. 




Figure 4: Energy dependence of the 
cross section of e + e _ — > 7r + 7r~7r°. Solid 
line shows the VDM prediction. 



Figure 5: Energy dependence of the 
cross section of e + e _ — > KsKl. Solid 
dots are from SND. 



• Cross section measurement for e + e — ► K$Kl 

The cross section of e + e~ — > K$Kl reaction was measured in 1982 by OLYA in 
Novosibirsk [38] and DM1 in Orsay|3£| in the energy region E cm = 1.06 — 1.40 GeV 
and E cm = 1.40 — 2.20 GeV respectively. It is desirable to re-measure this channel 
since the accuracy reached by both experiments is poor. So far, 1.8 pb _1 data has 
been analyzed by SND, utilizing K$ — > 7r°7r from e + e~ —>■ KsKl- e + e~ — ► um — > 
7r 7r°7 is the main background source. In addition, cosmic rays and beam associated 
background also contribute. The cross section measured by SND is illustrated in 
Fig. |5| with solid dots. 

• e + e _ — > U7T + 7T~ {uo — > 7T + 7r~7r°), t]7t + 7t~ (jj — > 7T + 7r~7r°, or 77) 

Figs. |5] and [7| plot the cross section measured by CMD-2 for um + 7i^ and t]7t + tt~, 
together with the measurement by DM2. The new measurement significantly re- 
duced the uncertainties though the systematic errors are still as high as 15%. 
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Figure 6: Cross sections for the reactions Figure 7: Cross section for e + e 
e + e~ — > uj7i + 7i^(uj — > 7T + 7r _ 7r°). ?77r + 7r~ (77 — ► 7r + 7r~7r°, or 77). 

• e + e~ — > 7r + 7r _ 

The cross section of the process e + e~ — > tt + tt~ is given by 



°=^0l\FM\\ (6) 

where F^s) and ^ are the pion form factor at the cm energy yfs and the velocity of 
the pion respectively. A precision measurement of the pion form factor is necessary 
to determine the i?-values via an exclusive method. The relative uncertainty contri- 



bution to is dominated by the e + e — > tt + -k channel with y/s < 2 GeV [|T^, |41 



The famous E821 experiment at BNL[46] has measured to a precision of ~ 5 ppm 



and will further improve the accuracy to about 1 ppm. In order to compare a mea- 
surement with such a high accuracy with theory, the uncertainty in R should be 
below 0.5% in this energy region. Therefore a new measurement of the pion form 
factor with smaller uncertainty is important for the interpretation of the E821 mea- 
surement. 

The pion form factor was measured by the OLYA and CMD groups at VEPP-2M 
about twenty years ago |47| . Twenty- four points from 360 to 820 MeV were studied 



by CMD with a systematic uncertainty of about 2%. The OLYA measurement 
scanned from 640 to 1400 MeV with small steps, giving a systematic uncertainty 
from 4% at the p-meson peak to 15% at 1400 MeV. 

The pion form factor is one of the major experiments planned at CMD-2. A 
total of 128 energy points were scanned in the whole VEPP-2M energy region (0.36- 
1.38 GeV) in six runs performed from 1994 to 1998 The discussion here is 

based on data taken from the first 3 runs with 43 energy points ranging from 0.61- 
0.96 GeV. The small energy scan step, 0.01 GeV, in this energy region allows the 
calculation of the hadronic contribution in a model-independent way. In order to 
investigate the a;- meson parameters and the p — u interference, the energy steps were 
2-6 MeV in the energy region near the tu-meson. The beam energy was measured 
with the resonance depolarization technique for almost all the energy points, which 
significantly reduced the systematic error arising from the energy uncertainty. The 
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charged trigger makes use of the information from the drift chamber and the Z- 
chamber and requires at least one track. There was an additional trigger criteria 
for the energy points between 0.81 and 0.96 GeV, which asks for the total energy 
deposited in the calorimeter to be greater than 20-30 MeV. The neutral trigger, 
served for monitoring the trigger efficiency, is based on the information only from 
the calorimeter. 

The background is mainly from cosmic muons. Bhabha and dimuon production 
are also background sources. The shape of the energy deposition was carefully 
studied for the event separation and selection. An event vertex cut was applied to 
reject cosmic muons effectively. 

To account for the fact that the radiative correction for e + e _ — > 7r + 7r~ depends 
on the energy behavior of the cross section of e + e _ — > 7r + 7r~ itself, the radiative 
correction factor was calculated iteratively. The existing |.FV(s)| 2 data were used as 
the first iteration for the calculation. The values of (^(s)) 2 were found to be stable 
after three iterations. 

The corrections for the pion losses due to decays in flight and nuclear interaction, 
as well as the background from uj —>■ 3tt were done using Monte Carlo simulations. 

The total systematic uncertainty was estimated to be currently 1.5% and 1.7% 
for the energy region 0.78-0.784 GeV and 0.782-0.94 GeV respectively, and 1.4% for 
all the other points. 

Other than the leading contribution from p(770) and a; (782), the resonances 
p(1450) and p(1700) should be taken into account to describe the data for the deter- 
mination of the pion form factor. In addition, the model based on the Hidden Local 
Symmetry (HLS), which predicts a point-like coupling 77r + 7r~, can well describe the 
experimental data below 1 GeV. Both the Gounaris-Sakurai (GS) parameterization 
and the HLS parameterization approaches WA 051 were used to fit the form factor. 



Only the higher resonance p(1450) was taken into account in fitting the pion form 
factor in the relatively narrow energy region 0.61-0.96 GeV. Fig. § shows the fit of 
the pion form factor with CMD-2 94, 95 data according to GS and HLS models. 
Both theoretical curves are indistinguishable. 




Figure 8: pion form factor measured by CMD-2 
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The pion form factor fit of CMD-2 94, 95 data is summarized by Table |2|. PDG 
data is also shown for comparison. 



Table 2: The results of fit of the CMD-2 (94,95) pion form factor data by GS and 
HLS models 





GS model 


HLS model 


M p (MeV) 


775.28±0.61±0.20 


774.57±0.60±0.20 


T p (MeV) 


147.70±1.29±0.40 


147.65±1.38±0.20 


Br{uj -> 7r + 7r-),% 


1.31±0.23±0.02 


1.32±0.23±0.02 


/3(GS) 


-0.0849 ± 0.0053 ±0.0050 




a(HLS) 




2.381±0.016±0.016 


X 2 /n 


0.77 


0.78 



With the remaining data collected, CMD-2 hopes to reduce the systematic error 
presented here by a factor of two. To achieve this goal, a new approach for the 
calculation of the radiative correction must be developed. 



2.2 R scan with BESII at BEPC in Beijing 

With the upgraded machine and detector [M, |5U], O], the BES collaboration per- 
formed two scans to measure R in the energy region of 2-5 GeV in 1998 and 1999. 
The first run scanned 6 energy points covering the energy from 2.6 to 5 GeV in the 
continuum. Separated beam running at each energy point was carried out in order 



to subtract the beam associated background from the data |pl 



The second run scanned about 85 energy points in the energy region of 2-4.8 



GeV[|2j. To subtract beam associated background, separated beam running was 
done at 26 energy points and single beam running for both e~ and e + was done 
at 7 energy points distributed over the whole scanned energy region. Special runs 
were taken at the J/ip to determine the trigger efficiency. The J/ip and '0(25') 
resonances were scanned at the beginning and at the end of the R scan for the 
energy calibration. 

The R values from the BESII scan data are measured by observing the final 
hadronic events inclusively, i.e. the value of R is determined from the number of 
observed hadronic events (N^ s d ) by the relation 

*l-L.e had -e trg -(l + 5)> { } 

where Nf, g is the number of beam associated background events; J2i Nu, (I = e, fi, r) 
and iV 77 are the numbers of misidentified lepton-pairs from one-photon and two- 
photon processes events; L is the integrated luminosity; 5 is the radiative correction; 
€had is the detection efficiency for hadronic events and e trg represents the trigger 
efficiency. 
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The trigger efficiencies are measured by comparing the responses to different 
trigger requirements in special runs taken at the J/ip resonance. From the trigger 
measurements, the efficiencies for Bhabha, dimuon and hadronic events are deter- 
mined to be 99.96%, 99.33% and 99.76%, respectively. As a cross check, the trigger 
information from the 2.6 and 3.55 GeV data samples is used to provide an indepen- 
dent measurement of the trigger efficiencies. This measurement is consistent with 
the efficiencies determined from the J/ip data. The errors in the trigger efficiencies 
for Bhabha and hadronic events are less than 0.5%. 

The task of the hadronic event selection is to identify one photon multi-hadron 
production from all other possible contamination mechanisms. The event selection 
makes full use of all the information from each sub-detector of BESII, namely, the 
vertex position, the measured charged-particle momentum and the energy loss due 
to ionization, the related time of flight, the associated pulse height and pulse shape 
of the electromagnetic calorimeter, and the hits in \x counter. 

The backgrounds involved in the measurement are mainly from cosmic rays, lep- 
ton pair production (e + e~, e + e~, r + r~), two-photon processes, and beam 

associated processes. The cosmic rays and part of the lepton pair production events 
are directly removed by the event selection. The remaining background from lep- 
ton pair production and two-photon processes is then subtracted out statistically 
according to a Monte Carlo simulation. 

The beam associated background sources are complicated. They may mainly 
come from beam-gas, beam-wall interaction, synchrotron radiation, and lost beam 
particles. The salient features of the beam associated background are that their 
tracks are very much along the beam pipe direction, the energy deposited in BSC is 
small, and most of the tracks are protons. 

Separated-beam runs were performed for the subtraction of beam associated 
background. Most of the beam associated background events are rejected by ver- 
tex and energy cuts. Applying the same hadronic events selection criteria to the 
separated-beam data, one can obtain the number of separated-beam events N sep 
surviving these criteria. The number of beam associated events N^ g in the corre- 
sponding hadronic events sample is given by Nb g = f ■ N sep , where / is the ratio of 
the product of the pressure at the collision region times the integrated beam currents 
for colliding beam runs and that for the separated beam runs. To subtract beam 
associated background in this way, the variation of the pressure in the collision re- 
gion and the beam current must be recorded for both colliding and separated-beam 
runs at each energy to be measured. 

JETSET7.4 is used as the hadronic event generator to determine the detection 
efficiency for hadronic events. Parameters in the generator are tuned using a 4 x 
10 4 hadronic event sample collected near 3.55 GeV for the tau mass measurement 



done by the BES collaboration [j5"5[| . The parameters of the generator are adjusted 
to reproduce distributions of kinematic variables such as multiplicity, sphericity, 
transverse momentum, etc. 

The parameters have also been obtained using the 2.6 GeV data (k 5 x 10 3 
events). The difference between the two parameter sets and between the data and the 
Monte Carlo data based on these parameter sets is used to determine a systematic 
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error of 1.9-3.2% in the hadronic efficiency. 

The Monte Carlo simulation packet JETSET was not designed to fully describe 
few body states produced by e + e~ annihilation in the few GeV energy region, though 
the event shapes are consistent with that from the Monte Carlo simulation with 
tuned parameters at 3.5 GeV. A great effort has been made by the Lund group 
and BES collaboration to develop the formalism using the basic Lund Model area 
law directly for the Monte Carlo simulation, which is expected to describe the data 
better |27|. 



Radiative corrections determined using four different schemes |p6| , |57| , |58| , ^£ 
agreed with each other to within 1% below charm threshold. Above charm threshold, 
where resonances are important, the agreement is within 1 ~ 3%. The major un- 
certainties common to all models are due to errors in previously measured i?-values 
and in the choice of values for the resonance parameters. For the measurements 
reported here, we use the formalism of ref. |58| and include the differences with the 
other schemes in the systematic error of 2.2-4.1%. 

The R values obtained at the six energy points scanned in 1998 are shown in 
Table 3 and graphically displayed in Fig. |9] with solid dots. The largest systematic 
error is due to the hadronic event selection and is determined to be 3.8-6.0% by 
varying the selection criteria. The systematic errors on the measurements below 4.0 
GeV are similar and are a measure of the amount of error common to all points. 
The BES collaboration has also performed the analysis including only events with 
greater than two charged tracks; although the statistics are smaller, the results 
obtained agree well with the results shown here. 



Table 3: -R-values measured with BESII from 1998 run. 



E cm (GeV) 


2.6 


3.2 


3.4 


3.55 


4.6 


5.00 


R 


2.64 


2.21 


2.38 


2.23 


3.58 


3.47 


Stat, error 


0.05 


0.07 


0.07 


0.06 


0.20 


0.32 


Sys. error 


0.19 


0.13 


0.16 


0.16 


0.29 


0.29 
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Figure 9: Plot of R values vs E c 
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The R values for E cm below 4 GeV are in good agreement with results from 772 
[0] and Pluto |9j but are below those from Mark I ||. Above 4 GeV, our values are 
consistent with previous measurements. 

Preliminary ^-values at 2.4, 2.5, 2.6, 2.7, 2.8, 2.9, 3.0, 4.6 and 4.8 GeV are plotted 
with solid squares in Fig. [l(| The preliminary errors, which add the statistical 
and systematical errors in quadrature, are all conservatively assigned to be 10%. 
However, it is believed that these errors can be decreased to be comparable to the 
error bars of the solid dots, i.e. ~ 7% for the energy points below 3.6 GeV and 
~ 10% for energies above 4.5 GeV. The first scan repeated 3.4 GeV and the second 
scan repeated 2.6 and 4.6 GeV data points measured in the first scan. In all case, 
i?-values obtained are consistent with each other at the same energy. 
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Figure 10: Energy points scanned with BESII. R values shown with triangles are 
from the online data estimation, of which the detection efficiency, the raditive cor- 
rection and the background have not yet been taken into account. The solid sqaures 
represent the preliminary i?-values obtained from BESII data, of which the total 
errors are all conservatively assigned to be 10%. The final i?-values of the 6 points 
listed in table 3 are also plotted with solid dots for comparison. 



3 Prospects and Concluding Remarks 

SND and CMD-2 at VEPP-2M have significantly improved the measurements of 
the hadron production cross section via e + e~ collisions for some of the important 
exclusive channels in the energy region of 0.36-1.38 GeV. Further improvement with 
the analysis of the existing data is forthcoming. A major advance would be possible 
if the energy region could be extended to 2 GeV, which would link up to the lowest 
energy of BEPC. 

CMD-2 and SND at VEPP-2M are planning to scan from threshold to 1.4 GeV 
in 1999-2000. A R scan between 2 to 10 GeV with KEDR at VEPP-4 is proposed. A 
scan covering such a wide energy region with the same machine and detector would 
be very important if the measurement could be performed with a ~ 1% precision. 
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The R scan performed with BESII at BEPC in Beijing can significantly reduce 
the uncertainties in R in the energy region 2-5 GeV. The i?-values from the first 
run data have already reduced the uncertainties in R from 15-20% to 7%. BESII 
at BEPC in Beijing is analyzing the second run R scan data. The preliminary R 
values in the whole energy region of 2-5 GeV are expected to be shown in the Spring 
of 2000, and the final results will be presented in the summer of 2000. 

The new R ratio results in e + e~ annihilations presented from Novosibirsk and 
Beijing have had a great impact on the value of ck(M§). Using these new (albeit, 
preliminary) results, A.D. Martin et. al. |Kj re-evaluate a(M§) and find a(M|) _1 = 



128.973 ± 0.035 or 128.934 ±0.040, according to whether inclusive or exclusive cross 
sections are used. The uncertainties here are already decreased more than half 
if we compare with the previous value of ct(Mf ) _1 = 128.89 ± 0.09 evalued by 



using the old experimental i?-values [22]. With the final results from the Beijing 



inclusive measurement in the whole 2-5 GeV region and the more precise results 
from Novosibirsk, a(M§) and a M will be even more precisely determined from the 
experimental data. 

A dedicated energy scan, aimed at a 1% precision direct measurement below 
1.4 GeV, planned by KLOE at DAPNE is not possible in the short term since 
the DAFNE machine is tuned for the resonance. However a machine upgrade is 
foreseen which will hopefully permit such an energy scan around 2004. Another 
method to measure the hadronic cross section is to measure events with Initial 
State Radiation (ISR). In this case one of the electrons or positrons of the beam 
radiates a hard photon and the cm energy of the hadronic system in the final state, 
mostly pions coming from the p-resonance, is lowered. KLOE has already started 



the analysis of those events [61 



Being one of the most fundamental parameters in particle physics, the i?-value 
plays an important role in the development of the theory of particle physics and in 
testing the Standard Model. Experimental efforts to precisely measure i?-values at 
low energies are crucial for the future electroweak precision physics. The measure- 
ments are not only important for the evaluation of a(M§) and for the interpretation 
of a but also necessary for the understanding of the hadron production mechanism 
via e + e~ annihilation. 

A real breakthrough in electroweak theory physics with regard to the .R-values 
at low energy would be possible only by measuring a(e + e~ — > hadrons) at ~ 1% 
accuracy. Such a level of precision requires significant improvement to both machine 
and detector, and needs better theoretical calculation of the radiative correction and 
the event generator for hadron production. 

Once the R- values has been measured with a precision of 1% in the energies 
covered by VEPP-2M and DAPNE, the central question will be how to further 
decrease the uncertainties of i?-values measured by BESII in the energy region of 
2-5 GeV, particularly from 2-3.7 GeV. Such a measurement would be extremely 
important for the interpretation of the experiment carrying out by E821 at BNL 
and for the precision determination of a(M§). This measurement would then be an 
important and attractive physics program for a r-c factory. 

I would like to thank the many people who have at some point or another helped 
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realize this presentation. In particular, my very special thanks are due Dr. Boris 
I. Khazin, Dr. Sergey Serednjakov for supplying me useful informations concerning 
the results from CMD-2 and SND. I'm also gratful to Dr. W. Dunwoodie, Dr. A. 
F. Haris and Dr. D. Kong who gave me helpful comments and suggestions. 
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